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摘 要 : 基于 同 古 孜 洛克 和 乌 鲁 瓦 提 水 文 站 1960 一 2016 年 年 尺度 与 月 尺度 径流 数据 ,分 别 采 用 均值 标准 差 法 和 


模糊 集 对 分 析 法 对 和 田 河 2 大 支流 进行 丰 枯 类 别 评价 ， 


= 


在 此 基础 上 优选 模糊 集 对 法 和 Copula 函数 法 计算 2 条 源 
流 径流 丰 枯 组 合 概率 ,同时 分 析 2 支 流 间 互补 的 优 劣 性 。 结 果 表 明 : 模 糊 集 对 法 考虑 了 年 内 径流 贡献 度 使 评价 结果 
更 客观 ,玉龙 喀什 河 枯 水 年 占 比 最 高 ,应 提高 枯水期 水 资源 利用 , 喀 拉 喀什 河 平水 年 占 比 最 高 ,表明 喀 拉 喀什 河 常 
FE 来 水 稳定 ;统计 法 和 Copula 函数 得 出 的 丰 枯 概率 存在 差异 性 与 数据 使 用 、 划 分 方式 及 源流 区 域 海拔 有 关 ,2 种 方 
法 结果 均 表 现 为 异步 丰 村 概率 大 于 同步 丰 村 概率 ,两 支流 具有 较 好 互补 性 。 该 结果 能 更 清晰 的 描述 和 田 河 径流 丰 
村 状态 及 丰 枯 组 合 特性 ,可 为 和 田 河 流域 水 资源 调度 提供 参考 和 决策 依据 。 


关键 词 : 模糊 集 对 法 ; Copula 函数 ; 均值 标准 差 法 ; 丰 桂 组合; 和 田 河 ; 新 疆 


河流 资源 作为 维系 自然 稳定 ,保障 生态 需 水 的 
重要 元 素 , 是 人 类 生产 活动 中 不 可 或 缺 的 资源 , 特 
别 是 在 西北 干旱 区 ,河川 径流 对 缓解 绿洲 生态 系统 
水 资源 短缺 的 矛盾 具有 重要 的 作用 。 近 年 来 , 受 全 
球 气 候 变 化 影响 ,和 田 河流 域 丰 枯 变化 规律 也 随 之 
改变 ,该 变化 直接 影响 了 和 田 河 流域 生态 与 社会 发 
展 需 水 的 变换 过 程 。 因 此 ,人 研究 环境 变化 下 和 田 河 
流域 变化 规律 ,对 干旱 区 水 资源 的 管理 提供 一 定 的 
参考 价值 。 

径流 丰 枯 划分 是 水 文 水 资源 计算 分 析 中 的 重 
要 内 容 , 常 用 方法 包括 水 文 频率 分 析 法 、 灰 色 分 类 、 
投影 寻 踪 法 .模糊 聚 类 法 、 集 对 分 析 法 等 4“。 上 述 方 
法 均 可 人 处理 水 文 分 类 问题 ,但 存在 忽略 径流 年 内 分 
配 , 求 解 复杂 及 具有 主观 判断 等 不 足 现象 “”。 虽 然 
集 对 分 析 法 弥补 了 上 述 不 足 ,但 未 考虑 边界 模糊 性 
及 年 内 径流 贡献 程度 问题 *”。 因 此 ,需要 引入 模糊 
理论 和 权重 分 析 构 建 模糊 集 对 方程 ,使 得 评价 结 
更 为 客观 准确 。 

丰 枯 变化 具有 一 定 的 规律 ,其 规律 性 表现 出 : 
在 相同 时 间 条 件 下 ,多 种 水 文 要 素 的 产生 丰 、 桂 情 
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况 的 不 同 , 即 丰 枯 组 合 概率 不 同 。 通 过 水 文 年 鉴 及 
文献 资料 2 统计 可 知 , 和 田 河 2 支流 (玉龙 喀什 河 与 
喀 拉 喀什 河 ) 出 现 的 丰 村 组 合 在 时 间 上 存在 重 蔷 部 
分 ,因此 ,本 文 针 对 变化 环境 下 的 喀 拉 喀什 河 、 玉 龙 
喀什 河 1960 一 2016 年 径流 量 ,综合 采用 模糊 理论 和 
权重 分 析 构 建 模糊 集 对 方程 ,对 2 支流 丰 枯 组 合 概 
率 进行 分 析 ,并 对 比 Copula 函数 构建 出 的 2 支流 间 
二 维 丰 村 组 合 概率 ,以 期 为 协调 和 田 河 流域 生产 、 
生活 、. 生 态 ( 三 生 ”) 用 水 ,缓解 水 资源 供给 矛盾 , 研 
究 径 流 丰 枯 变化 规律 ,评估 径流 丰 枯 类 别 可 为 流域 
内 洪 旱 灾害 的 防治 及 水 资源 高 效 利 用 提供 参考 
依据 。 


1 研究 区 概况 


和 田 河流 域 属 新 疆 西 南 区 域 ,位 于 塔里木 盆地 南 
缘 和 昆仑 山北 翡 之 间 , 地 理 坐 标 为 77"25'~81?"43' 下 、 
34°52'~40°28'N ,流域 面积 约 48870 km 。 受 昆仑 山 
和 帕 米 尔 高 原 阻 碍 影响 ,印度 洋 暧 湿 气 流 难以 进入 
流域 内 ,使 得 流域 内 形成 干旱 少雨 ,气温 高 ,潜在 
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蒸发 大 等 大 陆 性 气候。 和 田 河 干流 由 玉龙 喀什 河 
(简称 玉河 ) 与 喀 拉 喀什 河 (简称 喀 河 ) 汇 聚 形成 (图 
1) , 流 经 塔克拉玛干 沙漠 最 终 汇 入 塔里木 河 成 为 塔 
里 木 河 “ 三 源 " 之 一 ”。 和 田 河 不 仅 承担 向 塔里木 
河 输 水 任务 ,同时 也 维系 着 和 田 河 周边 地 区 的 生态 
环境 。 随 着 南 疆 地 区 基础 设施 不 断 完 善 与 经 济 投 
入 ,和 田 河 流域 人 口 基数 与 耕地 面积 快速 上 升 , 生 
态 用 水 受 生 活 .生产 用 水 挤占 ,流域 内 的 生态 稳定 
性 遭 到 威胁 , 硅 不 及 时 合理 分 配水 资源 利用 ,将 造 
成 流域 内 土地 荒漠 化 .生态 失衡 等 问题 ”。 


图 1 和 田 河 水 系 示意 图 
Fig. 1 Hotan River basin 


2 数据 与 方法 


在 传统 分 析 中 , 常 以 均值 标准 差 法 评判 径流 丰 
村 类 别 ,该 方法 是 以 年 径流 数据 为 基础 进行 丰 枯 判 
别 , 忽 略 了 年 内 时 段 径流 量 对 年 径流 总 量 贡献 的 不 
同 ,同时 该 方法 未 考虑 丰 枯 间 的 模糊 性 问题 易 导 致 
评价 结果 偏颇 ,因此 将 模糊 概念 与 权重 计算 加 入 集 
对 分 析 中 ,采用 模糊 集 对 法 对 径流 丰 枯 等 级 进行 科 
学 划分 。 对 于 径流 丰 村 变化 规律 的 研究 以 统计 法 
分 析 较 多 ,而 此 方法 需要 大 量 数据 作为 支撑 ,以 保 
证 研究 规律 的 精准 性 。 采 用 模糊 集 对 法 判断 径流 
丰 枯 类 别 ,并 在 此 评价 结果 上 ,由 统计 法 计算 出 丰 
村 变化 概率 。 采 用 Copula 函数 构建 丰 枯 组 合 模型 ， 
用 以 对 比 统计 法 结果 ,以 期 探究 径流 丰 枯 变化 
规律 。 
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研究 数据 采用 同 古 孜 洛克 和 乌 鲁 瓦 提 水 文 站 
1960 一 2016 年 实测 逐 月 径流 资料 进行 整理 分 析 。 
2.1 均值 标准 差 法 

以 年 尺度 径流 数据 为 基础 ,通过 相关 公式 划分 
5 种 丰 村 类 别 区 间 , 即 :| 0, + ks) 、[ E+ ksp, + 
kas;)s LE + hi Z tkis) [E tkis T+ kas) [E+ 
hs,+%), 其 中 x 和 ;为 年 平均 流量 的 均值 标准 
差 ,系数 经 验 取 值 分 别 为 =-0.9, 有 ;=-0.3, 友 =0.3， 
k.=0.9""), 

2.2 模糊 集 对 法 

以 月 径流 量 为 基础 ,通过 构建 模糊 集 对 式 判 断 
玉河 与 喀 河 径流 的 丰 枯 等 级 。 首 先 ,以 均值 标准 差 
法 确定 5 种 等 级 的 区 间 范 围 ”。 由 于 等 级 标准 符合 
越 小 越 优 的 条 件 约束 ,因此 采用 偏 小 型 也 数 进行 模 
糊 化 并 建立 联系 度 公 式 , 如 公式 (1) 所 示 。 


HL- = 

1 x SS 

S, +S, -2x 2x, ; 725, +58 
so ss ees 

S, +8, = 2a; ; 2x; a A DaT P < Sits 
53-5; Ss =: 2 

S,+5S,- 2x, , 2x, ,-8,-S,_ StS, Sats; 

<x. < 一 一 一 


S-S, S4-S, 2 meee 2 
2S, — 2x, ; p aTh Sa S345, < is Si 
S,-S; j S47 S3 l 2 i 
J ,Xi j> Sa 


(1) 

式 中 :7 为 差异 不 确定 分 量 系数 , 按 经 验 取 值 法 ,可 

Wize L=-0.5,L=0,4=0.5, XTS AM J=-1, RA BUR 

值 范围 为 -1 < py gp SV. HA as 为 集合 4i(i=1， 

2,…,n) 与 集合 B,(m=1,2,3,4,5) 的 联系 度 ,Si~5, 
为 临界 值 ,x;j 为 样本 值 。 

采用 炉 权 法 判断 年 内 径流 量 贡 献 比率 ,并 将 

计算 结果 代入 集 对 方程 中 ,得 到 综合 联系 度 公式 

(2); 
=> Ama s, = > ee Ajb; t+ 
a ar (2) 
SA b, sh tS Ab ait > Aged 

SUF sabi bobs ce NRA IE at n AWE A hs 

数 ; A, 为 第 j 个 指标 对 应 的 权重 值 ;ys 为 集合 联 

为 避免 系数 选择 存在 主观 干扰 ,采用 置信 和 度 准 
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则 对 样本 所 属 等 级 进行 判断 。 如 公式 (3) 所 示 : 
h,=f,tfitertf,>a (3) 


式 中 : /= 六 Aio ,为 第 个 指标 对 应 联系 度 分 量 ,n 


为 样本 总 数 ; 为 样本 所 属 第 级 ;a 为 置信 和 度 , 取 
值 界定 为 [0.5,0.7] °°’. 
2.3 Copula 函数 联合 分 布 

Copula 冰 数 是 连接 多 个 变量 边缘 分 布 所 构成 的 
多 维 联合 模型 ”” 。 其 假设 单 变量 边缘 分 布 分 别 为 
F(x) 有 0), 则 对 应 存在 二 维 联合 分 布 F(x,y) ,Copu- 
la 函数 则 可 以 表示 为 : 

F(x, y)=C[F (x), F,0)|= Clu, ») (4) 

式 中 ;uw 为 单 变量 的 边缘 分 布 。 

根据 函数 表达 式 进行 建 模 ,其 步骤 主要 为 : (1) 
通过 AIC(Akaike information criterion ) 检 验 和 均 方 根 
误差 (RMSE) 检 验 对 P-II 分布. 广 义 极 值 分 布 
(GEV) ,指数 分 布 (EXP) 和 正 态 分 布 (Normal) 进 行 
拟 合 优选 ,.(2) 通过 Kendall 相关 性 检验 计算 Copula 
函数 参数 。(3) 以 AIC. RMSE 检验 Gumbel, Frank, 
Gaussian 和 Clayton 函数 与 二 维 累 计 经 验 频率 拟 合 
效果 。(4) 以 累积 概率 P=62.5% 和 Pi=37.5% 为 让 


水 .平水 与 枯 水 的 分 割 频率 *!。 
3 实例 分 析 


3.1 径流 丰 枯 分 类 

3.1.1 基于 均值 标准 差 法 径流 丰 枯 分 类 ”根据 均值 
标准 差 法 对 径流 进行 划分 ( 表 1)。 

3.1.2 基于 模糊 集 对 法 径流 丰 枯 分 类 ”建立 玉河 、 
喀 河 1960 一 2016 年 月 径流 样本 xij(i=1, 2,*……,57; 
JH, 2,…, 12) ,并 确定 玉河 和 喀 河 不 同 指标 数 的 临 
界 值 , 表 2 为 不 同文 流 各 月 对 应 的 临界 值 。 

将 玉河 和 喀 河 不 同时 段 的 临界 值 和 x; 带 入 式 
(1) ,得 到 模糊 处 理 后 的 联系 度 分 量 。 以 玉河 1962 
年 1 月 径流 为 例 , 对 应 样本 为 xse=7.600, 介 于 Sis 
2oe<(Si+S)/2 ,选取 该 条 件 对 应 的 公式 计算 可 知 ， 
样本 xi 的 联系 度 分 量 包 隶属 枯 水 年 的 程度 为 
0.614 ,联系 度 分 量 w 属于 特 枯 水 年 的 程度 为 0.396， 
同 理 , 可 求 其 他 时 段 的 联系 度 分 量 。 

采用 炉 权 法 计算 玉河 、 喀 河 各 月 权重 ,可 知 玉 
河 权 重 值 A=(0.08245, 0.04707, 0.02942, 0.06947, 
0.09794, 0.13113, 0.10996, 0.14607, 0.06430, 
0.04725, 0.04512) ; 1% In] 4 H A Œ (E A =(0.03875, 


表 1 年 平均 径流 分 类 标准 


Tab.1 Classification standard of annual average runoff /(m*+s) 

类 见 特 村 村 F F EF 

玉河 <57.554 <67.883 <78.111 <88.390 > 88.390 

喀 河 <53.288 <65.980 <78.673 <91.365 291.365 

表 2 ERSA MEA 
Tab. 2 Classification standard of monthly runoff 
AO 玉河 喀 河 
Sı S2 Ss Ss Sı S2 Ss Ss 

1 6.526 8.026 9.525 11.025 8.068 9.249 10.429 11.609 
2 7.002 8.148 9.293 10.438 9.841 10.481 11.122 11.763 
3 7.687 8.625 9.561 10.499 10.382 13.317 16.252 19.186 
4 9.840 11.990 14.140 16.289 15.701 21.165 26.628 32.092 
5 22.703 28.603 34.502 40.402 36.532 45.996 55.460 64.924 
6 72.108 95.649 119.190 142.731 96.621 118.051 139.482 160.914 
7 192.605 254.435 316.266 378.097 163.288 210.746 258.205 305.663 
8 206.126 262.214 318.301 374.389 172.332 211.367 250.402 289.436 
9 47.396 64.683 81.970 99.256 50.609 68.568 86.527 104.860 
10 14.456 17.389 20.320 23.252 20.387 24.485 28.583 32.681 
11 9.229 10.775 12.282 13.806 13.313 16.126 18.939 21.753 
12 7.883 9.131 10.376 11.622 7.719 10.853 13.995 17.133 
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0.02170, 0.08951, 0.11829, 0.09184, 0.07064, 
0.10484, 1.33360, 0.12966, 0.05954, 0.06215, 
0.07943) 。 将 权重 值 联系 度 分 量 和 差异 不 确定 分 
量 系 数 代 入 公式 (2) 中 ,可 求 得 样本 隶属 某 评价 标 
准 的 可 能 性 。 样 本 wo 隶属 各 评价 等 级 结果 分 别 
为 :f=0.27756,f=0.48061 ,A=0.11233 ,fi=0.09542,, f= 
0.00349。 依 据 置信 度 准 则 筛选 样本 隶属 等 级 , 取 w= 
0.5, 可知 f+f>a, 由 此 判断 玉河 1962 年 为 村 水 年 , 即 
se 属于 开 类 标准 。 模 糊 集 对 法 分 类 结果 见 表 3。 
JAR 3 可知, 在 均值 标准 差 法 判断 下 ,玉河 近 
57a 内 共 出 现 15 次 丰 水 年 ,20 次 村 水 年 , 喀 河 近 57 a 
内 共 出 现 18 次 丰 水 年 ,26 次 枯 水 年 ;模糊 集 对 分 类 
结果 表明 :玉河 丰 水 年 有 10 次 , 枯 水 年 有 26 次 , 平 
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水 年 13 次 , 喀 河 丰 水 年 出 现 9 次 ,村 水 年 18 次 及 30 
次 平水 年 。 玉 河 、 喀 河 径 流 丰 枯 均 表现 出 , 枯 水 年 
次 数 大 于 丰 水 年 次 数 ,符合 干旱 区 水 资源 较 少 现 
象 。 这 表明 应 重点 关注 玉河 流域 枯 水 年 的 水 资源 
利用 ,而 喀 河 来 水 较为 平稳 。 
3.1.3 丰 枯 差异 分 析 对 比 2 种 丰 枯 评价 方法 结果 
可 知 , 玉 河 1960 年 .1971 年 .1973 年 .1994 年 和 2004 
年 . 喀 河 1961 年 .1966 年 和 2014 年 为 差异 度 最 大 年 
份 ,将 各 时 段 的 临界 值 绘 图 连 线 作 为 临界 线 ,以 玉 
河 1960 年 .1973 年 .1994 年 .2004 年 为 差异 代表 年 ， 
分 析 代 表 年 的 径流 过 程 (图 2)。 

由 图 2 分 析 1960 年 径流 过 程 ,1 一 6 月 和 10 一 12 
月 的 径流 均 在 S$, 以 下 ,径流 属于 “村 ”和 “ 特 枯 ”2 种 


RI 径流 丰 枯 分 类 结果 


Tab.3 Results of runoff classification 


年 从 模糊 集 对 法 均值 标准 差 法 年 从 模糊 集 对 法 均值 标准 差 法 

玉河 喀 河 玉河 喀 河 玉河 喀 河 玉河 喀 河 
1960 特 枯 村 丰 平 1989 村 村 特 枯 Ai 
1961 Æ 平 FE 特 丰 1990 平 平 平 丰 
1962 tii 特 村 枯 特 枯 1991 村 平 村 村 
1963 村 村 tii 枯 1992 特 枯 村 特 枯 枯 
1964 Hii Hii 村 村 1993 特 枯 枯 特 枯 特 村 
1965 特 枯 特 村 特 村 特 枯 1994 平 平 丰 
1966 £ thi E Æ 1995 Æ 村 村 ii 
1967 半 Æ Æ Æ 1996 £ £ 平 村 
1968 村 平 平 tii 1997 村 平 村 É 
1969 Ai 机 机 of 1998 F F F 村 
1970 平 平 平 平 1999 ar 村 aE 村 
1971 平 J 特 丰 平 2000 E Æ =e EF 
1972 ii 村 村 Ki 2001 平 平 平 Y 
1973 村 平 丰 平 2002 F F F 村 
1974 fii fii fii fii 2003 P F T P 
1975 Hii 枯 枯 Hih 2004 平 平 平 E 
1976 ii 村 村 Hih 2005 Æ £ £ £ 
1977 Æ 平 平 丰 2006 丰 丰 丰 丰 
1978 平 丰 特 丰 特 丰 2007 平 平 af F 
1979 tii 平 村 村 2008 平 平 Æ ii 
1980 村 下 枯 Hih 2009 Hii 平 村 村 
1981 Æ F 平 平 2010 FE $ 特 特 丰 
1982 Ai 平 平 平 2011 FE 丰 特 丰 
1983 平 Æ hii Æ 2012 FE 特 丰 特 丰 丰 
1984 枯 we we I 2013 FE 平 特 丰 we 
1985 Hii 平 平 Hii 2014 平 特 丰 平 村 
1986 特 枯 村 平 村 2015 平 特 丰 丰 特 
1987 村 村 特 村 村 2016 FE FE 特 丰 丰 
1988 村 平 平 = 
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图 2 代表 年 月 径流 过 程 分 布 
Fig. 2 Represents the distributions of annual and monthly 


runoff process 


状态 ,7 月 .9 月 径流 在 S$;、5: 之 间 , 分 类 结果 属于 “ 平 ” 
与 “ 丰 ” 之 间 , 由 于 径流 更 接近 5;, 属 于 “ 平 ” 的 概率 更 
大 ,因此 7 月 、9 月 分 类 等 级 为 “ 平 ”,8 月 径流 大 于 临 
界 值 S$, 属于 “ 特 丰 ”, 因 此 判定 玉河 1960 年 为 “ 特 
村 ”水 年 。 同 理 可 分 析 其 余 差 异 年 份 。 
差异 原因 在 于 模糊 集 对 法 考虑 到 年 内 径流 对 

全 年 径流 贡献 程度 不 同 , 而 玉河 与 喀 河 为 季节 性 河 
流下 ,具有 汛期 径流 量 大 , 非 汛 期 径流 量 少 的 特点 ， 
因此 在 赋 权 后 的 模糊 集 对 评价 结果 与 均值 标准 法 
评价 结果 存在 差异 。 而 均值 标准 差 法 视 各 时 有 段 径 
流 比 重 一 致 ,从 而 扩大 了 年 内 时 段 中 极端 情况 , 导 
致 评价 结果 偏 大 。 通 过 计算 可 知 ,玉河 1960 年 汛期 
径流 权重 为 51.6% , 非 汛期 径流 权重 为 49.495 ,两 者 
相差 不 大 ,这 表明 ,径流 丰 枯 等 级 判定 与 径流 量 
小 和 年 内 径流 的 权重 均 相 关 。 模 糊 集 对 法 能 综合 
考虑 时 程 分 配 , 能 客观 准确 的 表达 年 径流 分 类 
结果 。 
3.2 径流 丰 枯 组 合 分 析 
3.2.1 基于 模糊 集 对 分 析 结 果 的 丰 枯 组 合 概率 分 析 

上 述 丰 枯 评 价 方法 均 可 分 析 径 流 丰 枯 状 态 , 因 模 
糊 集 对 评价 法 能 全 面 .客观 反应 出 径流 真实 状态 ， 
以 模糊 集 对 评价 结果 为 基础 ,基于 玉河 、 喀 河 的 月 


概率 。 

为 方便 计算 ,将 径流 分 为 让、 平 . 村 3 类 进行 组 
合 概 率 分 析 ,将 模糊 集 对 法 分 类 结果 中 的 “ 特 丰 ” 归 
入 “ 丰 ”,“ 特 枯 ” 归 于 “村 ”, 则 2 河 共 形成 9 种 组 合 状 
态 ( 表 4)。 

玉河 和 喀 河 同步 丰 枯 概率 为 33.833% , 同 丰 的 概 
率 最 小 为 2.77%, 同 平 . 同 村 的 概率 分 别 为 17.54% 和 
15.51%。 异 步 丰 村 组 合 概率 为 64.17%, 大 于 同步 丰 
枯 组 合 概率 ,初步 可 断定 ,利于 河流 间 的 丰 枯 互 
补 。 玉 丰 喀 枯 的 组 合 概率 为 5.54% , 玉 枯 喀 丰 的 组 
合 概率 为 7.75% ,表明 玉河 径流 与 喀 河 径流 最 佳 丰 
村 互补 概率 较 低 。 

3.2.2 基于 Copula 函数 丰 枯 组 合 概率 分 析 Copula 
函数 是 以 玉河 ` 喀 河 年 径流 量 为 基础 ,将 其 序列 的 经 
验 频率 分 别 与 4 种 单 变 量 边缘 分 布 函数 拟 合 ( 表 5)。 

从 表 5 的 结果 可 知 ,在 AIC 检验 与 RMSE 检验 
下 ,玉河 、 喀 河 服从 P-II 分 布 检验 值 为 最 小 值 , 因 此 
以 P-II 分 布 作 为 2 径流 序列 的 最 优 边 缘分 布 构建 联 
合 分 布 。 

拟 合 二 维 联合 频率 和 经 验 频率 ,并 通过 AIC 和 
RMSE 对 拟 合 结果 进行 检验 。 在 4 种 Copula 函数 类 
型 中 ,Frank-Copula 函数 拟 合 效 果 最 好 ,AIC RMSE 

仿 验 结果 分 别 为 -429.895 ,0.022 。 

通过 Kendall 相关 计算 得 出 系数 r=0.6295 , 求 得 
Copula 参 数值 0=0.8354。 

根据 丰 枯 划分 频率 ,求解 出 Copula 函数 联合 丰 
村 组 合 概率 ( 表 6)。 

由 表 6 可 知 , 径 流 同 步 丰 枯 概 率 总 计 达 到 
48.3% , 其 中 同 丰 、 同 平 、 同 枯 的 概率 为 分 别 为 
8.41% 18.48% 、21.41% , 同 枯 组 合 概率 为 9 种 组 合 
状况 中 最 高 值 ,异步 丰 枯 概率 为 531.7% ,其 中 , 丰 - 枯 
最 有 利 互补 组 合 的 组 合 概率 为 16.66%。2 文 流 整体 
组 合 概率 表现 出 异步 组 合 概率 总 和 大 于 同步 丰 枯 
组 合 概率 ,说 明 2 支 流 具 有 良好 的 互补 性 , 且 该 结 
与 统计 法 得 到 的 组 合 概率 表现 一 致 ,说 明 Copula PK 


表 4 基于 统计 结果 下 的 丰 枯 组 合 概率 


Tab.4 Probability of synchronous-asynchronous encounters based on statistical results 


丰 桔 同步 概率 /% 


河流 


丰 村 异步 概率 /% 


同 丰 同 平 同村 合计 EER 


E EFR EPRE EPER RANE ERE 。 合计 


玉河 - 喀 河 2.77 17.54 15.51 35.83 9.23 


5.54 5.26 10.52 7.75 25.85 64.17 
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表 5 变量 拟 合 检 验 
Tab. 5 Fitness test results 
分 布 状态 河流 站 点 检验 方法 ie 

P- GEV EXP Normal 
单 变 量 分 布 玉河 AIC -405.534 -356.028 -392.461 -356.748 
RMSE 0.027 0.032 0.028 0.032 
喀 河 AIC -410.029 -332.626 -367.053 =332.017 
RMSE 0.026 0.034 0.03 0.035 

联合 变量 分 布 玉河 - 喀 河 Frank Gumbel Gaussian Clayton 
AIC —429.895 —418.378 —408.106 -397.119 

RMSE 0.022 0.024 0.027 0.03 


#6 基于 Copula 函数 的 丰 枯 组 合 概率 


Tab.6 Probability of synchronous-asynchronous encounters of the Copula joint distributions 


河流 丰 枯 同步 概率 /% 丰 枯 异步 概率 /% 
同 丰 同 平 同村 合计 玉 丰 喀 平 EER EPRE FOR MAE EMR ”合计 
玉河 - 喀 河 8.41 18.48 21.41 48.3 5.76 8.33 5.76 11.76 8.33 11.76 51.7 
数 分 析 结 果 合 理 。 径流 量 可 以 得 到 应 证 六 ,由 于 不 同 海拔 对 应 融雪 速 
or SEAS E, Nc wm EE LL ak SS Pah NY TE PB, 
4 讨论 冰雪 融化 时 间 较为 滞后 ,因此 造成 了 丰 枯 异步 概率 


对 比 模糊 集 对 结果 统计 概率 与 Copula 函数 联 
合 概率 ,两 者 结果 整体 相同 ,但 丰 枯 组 合 概率 有 所 
不 同 , 丰 枯 同步 最 大 相差 5.9% ,让 村 异步 最 大 相差 
14.09%。 因 此 ,分 析 其 造成 差异 原因 可 能 存在 :(1) 
Copula 函数 是 基于 年 尺度 数据 构建 联合 分 布 函 数 ， 
以 联合 概率 计算 出 丰 枯 组 合 概率 ;模糊 集 对 分 析 是 
以 月 尺度 数据 评价 出 年 径流 丰 枯 结果 ,利用 统计 方 
法 计算 出 丰 枯 组 合 频率 ,两 者 在 数据 使 用 上 有 所 不 
同 。(2) Copula 函数 计算 丰 村 组合 中 利用 37.5% 和 
62.5% 的 频率 划分 径流 丰 桂 存在 主观 判断 ,(3) 模糊 
集 对 结果 得 出 的 统计 频率 是 通过 近 57 a 的 径流 资料 
而 计算 出 了 组 合 频率 ,水 文 序列 长 度 有 限 ,不 能 代 
表 整 体 事件 的 概率 ,而 Copula 函数 能 从 已 知 的 二 维 
联合 频率 趋势 来 推 基 较 长 时 段 的 径流 丰 枯 组 合 的 
概率 。 

同时 借助 于 物理 成 因 分 析 可 知 , 玉 河 主要 产 流 
区 西 昆 仑 山 和 喀 河 产 流 区 喀 喇 昆仑 山 对 应 的 积 雪 
冰川 覆盖 面积 海拔 的 不 同 。 咯 喇 昆 仑 山 冰雪 覆盖 
面积 较 广 海拔 较 高 且 年 均 0 名 等 温 线 海拔 为 4200 
mm 左右 ; 西 昆 仑 山 冰 雪 履 善 面积 小 于 喀 喇 昆仑 山 ,其 
最 大 融雪 海拔 介 于 3800~4000 m Z Ù ,形成 有 效 径 
流 补 给 大 于 喀 喇 昆仑 山 冰 雪 补 给 ,这 点 通过 对 比 年 


较 大 的 特征 。 

为 提高 水 资源 利用 效率 与 经 济 发 展 , 乌 鲁 瓦 提 
水 库 作 为 重点 建设 项 目 于 2002 年 修建 在 喀 河 出 山 
口 处 ,其 控制 流域 面积 19983 km , 兴 利 库容 达 2.24x 
10m ,对 洪水 具有 显著 削 峰 作用 ,使 得 汛期 洪峰 流 
量 小 于 天 然 洪 峰 流 量 。 相 关 研 究 表明 , 乌 鲁 瓦 提 水 
库 的 修建 使 喀 河 年 内 径流 分 配 更 为 均匀 , 且 通 过 水 
库 提前 泄洪 排 沙 , 使 得 喀 河 与 玉河 的 同步 丰 顶 概率 
较 低 55。 本 文 将 模糊 集 对 评价 结果 与 基于 Copula 
函数 方法 计算 的 丰 枯 组 合 概率 进行 对 比 ,2 种 方法 
得 到 结果 均 可 以 反应 水 文 事件 概率 。 在 短期 水 文 
序列 研究 中 ,采用 统计 法 得 到 的 结果 更 易 反映 近年 
来 流域 丰 枯 的 组 合 状况 ,为 水 资源 合理 优化 配置 提 
供 参 考 依 据 。 而 从 水 利 工程 设计 角度 出 发 ,采用 
Copula 函数 联合 分 布 得 到 的 组 合 概率 能 更 好 为 水 利 
工程 设计 提供 支撑 依据 。 同 时 本 文 所 探讨 的 方法 
可 为 径流 丰 枯 评价 及 组 合 分 析 提 供 更 好 的 研究 思 
路 ,为 丰 枯 组 合 概率 研究 提供 理论 参考 。 


5 结论 


径流 丰 枯 变化 规律 是 流域 水 资源 规划 配置 . 管 
理 的 主要 依据 。 人 研究 径流 直 村 评价 及 组 合 分 析 , 可 
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习 
由 


有 效 的 反映 径流 丰 枯 变化 情况 ,从 而 为 流域 水 资源 
管理 提供 理论 依据 。 本 文 基于 均值 标准 差 法 及 模 
糊 集 对 分 析 方 法 ,并 结合 Copula 函数 ,探求 和 田 河 
流域 丰 枯 组 合 概率 ,其 结论 如 下 : 

(1) 在 径流 丰 枯 分 类 中 ,均值 标准 差 法 以 年 径 
流 数据 为 基础 ,未 考虑 径流 年 内 变化 ,导致 评价 结 
果 较 为 片面 。 而 相 较 均值 标准 差 法 ,模糊 集 对 分 析 
方法 考虑 了 年 内 径流 分 布 和 临界 模糊 问题 对 径流 
丰 枯 分 类 的 影响 ,并 以 水 文 站 各 年 份 的 月 径流 构建 
模糊 联系 度 方程 进行 丰 枯 分 类 ,使 评价 结果 客观 
合理 。 

(2) 由 模糊 集 对 法 可 知 , 由 于 玉龙 喀什 河 在 枯 
水 年 占 比 较 高 ,应 提高 枯水期 水 资源 利用 量 ; 同时 ， 
喀 拉 喀什 河 平水 年 的 占 比 较 高 ,充分 表现 出 喀 拉 喀 
什 河 常年 来 水 较为 稳定 。 

(3) 通过 丰 枯 组 合 分 析 ,2 河 流 异 步 丰 枯 组 合 概 
率 大 于 同步 丰 枯 组 合 概率 ,径流 互补 性 较 好 , 且 发 
生 特大 洪水 概率 较 小 ,此 研究 可 为 流域 水 资源 管理 
提供 坚实 的 基础 。 
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Evaluation and combination analysis of runoff in Hotan River, Xinjiang 


HUANG Xing', CHEN Fulong', ZHAO Qi’, HE Chaofei', LONG Aihua’’ 


(1. College of Water Conservancy & Architectural Engineering, Shihezi University, Shihezi 832000, Xinjiang, China; 
2. Hydrology Bureau of Xinjiang Uygur Autonomous Region, Urumqi 830000, Xinjiang, China; 3. State Key 
Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and 
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Abstract: Studies on the different contributions of the runoff of the Hotan River in different time periods and 
synchronous- asynchronous encounter probability can provide a reference for the balance of industrial water, 
potable water, and water ecological environment and the improvement of water resource utilization efficiency in 
the Hotan River Basin. Mean standard deviation and fuzzy set pair analysis methods are applied to evaluate the 
runoff wetness-dryness of the Yulong Kashi River and the Karakashi River. The possibility of inter-runoff water 
transfer and synchronous-asynchronous encounter probability is analyzed on the basis of fuzzy set pair method 
and copula function. The fuzzy set pair method reveals that the Yulong Kashi River had 10 and 26 wet and dry 
years, respectively. By comparison, the mean standard deviation method shows that the Yulong Kashi River had 
15 and 20 wet and dry years, respectively. The Karakash River had 9 wet years. Furthermore, 18 and 26 times in 
dry years were observed. In the fuzzy set pair method, the contribution of annual runoff and the influence of time 
history distribution are considered. Thus, evaluation results are more objective, and the difference in the 
probability of abundance and dryness obtained by the statistical method and the copula function is the same as the 
data scale is related to the way of dividing the abundance and dryness. The results of the two methods 
demonstrate that the probability of asynchronous abundance and dryness is greater than that of synchronous 
abundance and dryness, and the two tributaries have good complementarity. They can more clearly describe the 
high and low runoff status and the characteristics of high and low runoff in the Hotan River. They can also 
provide a basis for making decisions on water resource dispatch in the Hotan River Basin. 

Keywords: fuzzy set pair method; Copula function; mean standard deviation method; synchronous-asynchro- 


nous combined; Hotan River; Xinjiang 


